o 
o 

(N 



Oh! 



SPIN-06/16 
ITP-UU-06/19 



Quantum Gravity, or 
The Art of Building Spacetime^ 

J. Ainbj0rn'^''^ J. Jurkiewicz'^ , and R. LoW^ 



"• The Niels Bohr Institute, Copenhagen University 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark. 
rS ,' emaih ambjorn@nbi.dk 

< 

(N 



^ Institute of Physics, Jagellonian University, 

Reymonta 4, PL 30-059 Krakow, Poland. 

email: jurkicwi@thrisc.if.uj.cdu.pl 



^>. , "^ Institute for Theoretical Physics, Utrecht University, 

^ ■ Leuvenlaan 4, NL-3584 CE Utrecht, The Netherlands. 

email: j.ambjorn@phys.uu.nl, r.loll@phys.uu.nl 

O 

o 



Abstract 



The method of four-dimensional Causal Dynamical Triangulations provides a 
background-independent definition of the sum over geometries in quantum grav- 
/\ I ity, in the presence of a positive cosmological constant. We present the evidence 

c^ ' accumulated to date that a macroscopic four- dimensional world can emerge from 

this theory dynamically. Using computer simulations we observe in the Euclidean 
sector a universe whose scale factor exhibits the same dynamics as that of the 
simplest mini-superspace models in quantum cosmology, with the distinction that 
in the case of causal dynamical triangulations the effective action for the scale 
factor is not put in by hand but obtained by integrating out in the quantum theory 
the full set of dynamical degrees of freedom except for the scale factor itself. 



^Contribution to the book " Approaches to Quantum Gravity", edited by D. Oriti, to appear 
at Cambridge University Press. 
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1 Introduction 

What is more natural than constructing space from elementary geometric building 
blocks? It is not as easy as one might think, based on our intuition of playing 
with Lego blocks in three-dimensional space. Imagine the building blocks are 
(i-dimensional flat simplices all of whose side lengths are a, and let d > 2. The 
problem is that if we glue such blocks together carelessly we will with probability 
one create a space of no extension, in which it is possible to get from one vertex to 
any other in a few steps, moving along the one-dimensional edges of the simplicial 
manifold we have created. We can also say that the space has an extension which 
remains at the 'cut-off' scale a. Our intuition coming from playing with Lego 
blocks is misleading here because it presupposes that the building blocks are 
embedded geometrically faithfully in Euclidean IR^, which is not the case for the 
intrinsic geometric construction of a simplicial space. 

By contrast, let us now be more careful in our construction work by assigning 
to a simplicial space T - which we will interpret as a (Euclidean) spacetime - 
the weight e"'^'-^-', where S{T) denotes the Einstein action associated with the 
piecewise linear geometry uniquely defined by our construction^. As long as the 
(bare) gravitational coupling constant Gn is large, we have the same situation 
as before. However, upon lowering Gn we will eventually encounter a phase 
transition beyond which the geometry is no longer crumpled into a tiny ball, but 
maximally extended. Such a geometry is made out of effectively one-dimensional 
filaments^ which can branch out, and are therefore called branched polymers or 
trees [H |2] . The transition separating the two phases Hj is of first order, which 
implies that there is no smooth change between the two pathological types of 
minimally or maximally extended "universes" . 

In order for the sum over geometries to produce a quantum theory of gravity 
in which classical geometry is reproduced in a suitable limit, we therefore need a 
different principle for selecting the geometries to be included in this sum. Below 
we will introduce such a principle: our prescription will be to sum over a class of 
(Euclidean) geometries which are in one-to-one correspondence with Lorentzian, 
causal geometries. At the discretized level, where we use a specific set of building 
blocks and gluing rules to constructively define the path integral, we call these 
geometries causal dynamical triangulations (CDT) jHl El 13 E] ■ 

Before discussing CDT in more detail let us comment on the nature of the 
geometries contributing to the path integral. It is important to emphasize that in 
a quantum theory of gravity a given spacetime geometry as such has no immediate 



^There exists a natural, coordinate-independent definition of the Einstein action for piecewise 
linear geometries called the Regge action. 

■'The rf-dimensional building blocks are arranged such that {d — 1) "transverse" dimensions 
have a size of only a few lattice spacings. 



physical meaning. The situation is really the same as in ordinary quantum field 
theory or even quantum mechanics, where individual field configurations 0(x, t) 
or particle paths x{t) are not observable. Only certain expectation values related 
to the fields or paths can be observed in experiments. This does not mean there 
cannot exist limits in which it is appropriate to talk about a particular field 
configuration or the path of a particle in an approximate sense. In the case 
of our actual universe, down to the smallest distances that have been probed 
experimentally, it certainly does seem adequate to talk about a fixed classical 
spacetime geometry. Nevertheless, at sufficiently small distances it will no longer 
make sense to ask classical questions about spacetime, at least if we are to believe 
in the principles of conventional quantum theory. 

By way of illustration let us discuss the situation for the ordinary harmonic os- 
cillator (or the free particle) and consider the path integral from (a;i, ti) to {x2, ^2). 
Precisely for the harmonic oscillator (or the free particle) the decomposition 

x{t) = x,i{t) + y{t), y{t,) = y{t2) = 0, (1) 

leads to an exact factorization of the path integral, because the action satisfies 

S{x) = S{x,i) + S{y). (2) 

This implies that the classical path Xci{t) contributes to the path integral with the 
classical action, and y{t) with quantum fiuctuations independent of this classical 
part. Taking the classical trajectory to be macroscopic one obtains the picture of 
a macroscopic path dressed with small quantum fiuctuations; small because they 
are independent of the classical motion. An explicit Euclidean calculation yields 
the result 

f dt y\t)\ = ^-^ (cuTtanh-^ ujT - 1) (3) 

as a function of the oscillator frequency u and mass m. Let us now consider 
a situation where we have chosen the "system size", i.e. Xci{t), to be macro- 
scopic. According to 0, the quantum fiuctuations around this path can then be 
considered small since h is small. 

This is more or less the picture we envisage for our present-day universe in 
quantum gravity: the universe is of macroscopic size, governed by the classical 
equations of motion (the analogue of choosing "by hand" {xi,ti) and (x2,t2) to 
be macroscopic in the example above), and the small quantum fiuctuations are 
dictated by the gravitational coupling constant (times h/c"). 

A given configuration x{t) in the path integral for the quantum-mechanical 
particle is (with probability one) a continuous, nowhere differentiable path, which 
moreover is fractal with Hausdorff dimension two, as we know from the rigorous 
construction of the Wiener measure on the set of parametrized paths. In the 



case of quantum gravity we do not have a similar mathematically rigorously de- 
fined measure on the space of geometries, but it is natural to expect that if it 
exists, a typical geometry in the path integral will be continuous, but nowhere 
differentiable. By analogy, the piecewise linear geometries seem a good choice if 
we want to approximate the gravitational path integral by a set of geometries 
and subsequently take a limit where the approximation (the cut-off) is removed. 
Moreover, such simplicial manifolds possess a natural, geometric and coordinate- 
independent implementation of the Einstein-Hilbert action. With all local curva- 
ture degrees of freedom present (albeit in a discretized fashion), we also expect 
them to be suitably "dense" in the set of all continuous geometries. 

The spirit is very much that of the standard lattice formulation of quantum 
field theory where (fiat) spacetime is approximated by a hypercubic lattice. The 
ultraviolet cut-off in such field theories is given by the lattice spacing, i.e. the 
length of all one-dimensional lattice edges. We can in a similar and simple manner 
introduce a diffeomorphism-invariant cut-off in the sum over the piecewise linear 
geometries by restricting it to the building blocks mentioned earlier. A natural 
building block for a d-dimensional spacetime is a d-dimensional equilateral sim- 
plex with side-length a, and the path integral is approximated by performing the 
sum over all geometries (of fixed topology"^) which can be obtained by gluing such 
building blocks together, each geometry weighted appropriately (for example, by 
e~^ , where S is the Einstein-Hilbert action). Afterwards we take the limit a — > 0. 
For a particular choice of the bare, dimensionless coupling constants one may be 
able to obtain a continuum limit, and thus extract a continuum theory. For other 
values, if the sum exists at all (possibly after renormalization), one will merely 
obtain a sum which has no continuum interpretation. This situation is precisely 
the same as encountered in ordinary lattice field theory in fiat spacetime. 

As mentioned earlier it has up to now not been possible to define construc- 



^In classical General Relativity there is no motivation to consider spacetimes whose spatial 
topology changes in time, since their Lorentzian structure is necessarily singular. There is an 
interesting and long-standing discussion about whether one should include topology changes in 
a quantum theory of gravity. However, even in the case of two-dimensional Euclidean quantum 
gravity, where the classification of topology changes is simple, the summation over topologies 
has never been defined non-perturbatively in a satisfactory way, despite many attempts, in 
particular, in so-called non-critical string theory. (However, see Pj for how one may improve 
the convergence of the sum in two-dimensional Lorentzian quantum gravity by invoking not 
just the topological, but the causal, geometric structure of spacetime.) The situation becomes 
worse in higher dimensions. For instance, four-dimensional topologies are not classifiable, so 
what does it mean to sum over them in the path integral? The problem - even in dimension 
two - is that there are many more geometries of complicated topology than there are of simple 
topology, with the consequence that any sum over geometries will be (i) completely dominated 
by these complicated topologies, and (ii) plainly divergent in a way which (until now) has 
made it impossible to define the theory non-perturbatively in an unambiguous and physically 
satisfactory manner. In higher dimensions these problems are totally out of control. 



tively a Euclidean path integral for gravity in four dimensions by following the 
philosophy just outlined. One simply has not succeeded in identifying a contin- 
uum limit of the (unrestricted) sum over Euclidean building blocks. Among the 
reasons that have been advanced to explain this failure, it is clear that the entropy 
of the various geometries plays an important role. We have already pointed out 
that the crumpled geometries of no extension dominate the space of all continu- 
ous geometries whenever the dimension of spacetime is larger than two. There is 
nothing wrong with this a priori; the path integral of any quantum field theory 
is dominated completely by wild UV-field fluctuations. However, in the case of 
renormalizable quantum field theories there exists a well-defined limiting proce- 
dure which allows one to extract "continuum" physics by fine-tuning the bare 
coupling constants of the theory. An analogous procedure in Euclidean quantum 
gravity still has not been found, and adding (bosonic) matter does not improve 
the situation. Instead, note that the Einstein-Hilbert action has a unique feature, 
namely, it is unbounded from below. The transition between the crumpled and 
the branched-polymer geometries can be seen as a transition from a phase where 
the entropy of configurations dominates over the action to a phase where the 
unboundedness of the Euclidean action becomes dominant^. The impossibility 
of finding a continuum limit may be seen as the impossibility of balancing the 
entropy of configurations against the action. We need another guiding princi- 
ple for selecting Euclidean geometries in the path integral in order to obtain a 
continuum limit, and it is such a principle we turn to next. 

2 Defining CDT 

It has been suggested that the signature of spacetime may be explained from a 
dynamical principle ^T]. Being somewhat less ambitious, we will assume it has 
Lorentzian signature and accordingly change our perspective from the Euclidean 
formulation of the path integral discussed in the previous section to a Lorentzian 
formulation, motivated by the uncontroversial fact that our universe has three 
space and one time dimension. A specific rotation to Euclidean signature intro- 
duced below will be needed in our set-up as a merely technical tool to perform 
certain sums over geometries. Unlike in fiat spacetime there are no general theo- 



^Although the action is not unbounded below in the regularized theory, this feature of the 
continuum action nevertheless manifests itself in the limit as the (discretized) volume of space- 
time is increased, eventually leading to the above-mentioned phase transition at a particular 
value of the bare gravitational coupling constant. Remarkably, a related phenomenon occurs 
in bosonic string theory. If the world-sheet theory is regularized non-perturbatively in terms 
of triangulations (with each two-dimensional world-sheet glued from fundamental simplicial 
building blocks), the tachyonic sickness of the theory manifests itself in the form of surfaces 
degenerating into branched polymers |10|. 



rems which would allow us to relate the Euclidean and Lorentzian quantum field 
theories when dealing with quantum gravity. 

Consider now a connected space-like hypersurface in spacetime. Any clas- 
sical evolution in general relativity will leave the topology of this hypersurface 
unchanged, since otherwise spacetime would contain regions where the metric is 
degenerate. However, as long as we do not have a consistent theory of quan- 
tum gravity we do not know whether such degenerate configurations should be 
included in the path integral. We have already argued that the inclusion of ar- 
bitrary spacetime topologies leads to a path integral that has little chance of 
making sense. One might still consider a situation where the overall topology of 
spacetime is fixed, but where one allows "baby universes" to branch off from the 
main universe, without permitting them to rejoin it and thus form "handles". 
Apart from being a rather artificial constraint on geometry, such a construction 
is unlikely to be compatible with unitarity. We will in the following take a con- 
servative point of view and only sum over geometries (with Lorentzian signature) 
which permit a foliation in (proper) time and are causally well-behaved in the 
sense that no topology changes are allowed as a function of time. In the context 
of a formal continuum path integral for gravity, similar ideas have earlier been 
advanced in [T2]. 

Of the diffeomorphism-invariant quantities one can consider in the quantum 
theory, we have chosen a particular proper-time propagator, which can be defined 
constructively in a transparent way. We are thus interested in defining the path 

integral 

r9{T) 
G{g{Ql g{T)-T) = Vg e^'^^^ (4) 

■^9(0) 

over Lorentzian geometries on a manifold Ai with topology S x [0, 1], where S 
is a compact, connected three-dimensional manifold. The geometries included in 
the path integral will be such that the induced boundary three-geometries g{0) 
and g{T) are space-like and separated by a time-like geodesic distance T, with T 
an external (diffeomorphism-invariant) parameter. 

We now turn to the constructive definition of this object in terms of build- 
ing blocks. The discretized analogue of an infinitesimal proper-time "sandwich" 
in the continuum will be a finite sandwich of thickness At = 1 (measured in 
"building block units" a) of topology S x [0, 1] consisting of a single layer of four- 
simplices. This layer has two spacelike boundaries, corresponding to two slices 
of constant (integer) "proper time" t which are one unit apart. They form two 
three-dimensional piecewise fiat manifolds of topology S and consist of purely 
spacelike tetrahedra. By construction, the sandwich interior contains no vertices, 
so that any one of the four-simplices shares k of its vertices with the initial spa- 
tial slice and 5 — k oi them with the final spatial slice, where 1 < k < 4. To 
obtain extended spacetimes, one glues together sandwiches pairwise along their 
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(4,1) (3,2) 

Figure 1: The two fundamental building blocks of causal dynamically triangulated 
gravity. The flat four-simplex of type (4,1) on the left has four of its vertices at time 
t and one at time t+1, and analogously for the (3,2)-simplex on the right. The "gap" 
between two consecutive spatial slices of constant integer time is filled by copies of 
these simplicial building blocks and their time-reversed counterparts, the (1,4)- and 
the (2,3)-simplices. 

matching three-dimensional boundary geometries. We choose each four-simplex 
to have time-like links of length-squared af and space-like links of length-squared 
a^, with all of the latter located in spatial slices of constant integer-t. 

Each spatial tetrahedron at time t is therefore shared by two four-simplices 
(said to be of type (1,4) and (4,1)) whose fifth vertex lies in the neighbouring slice 
of constant time t — 1 and t + 1 respectively. In addition we need four-simplices 
of type (2,3) and (3,2) which share one link and one triangle with two adjacent 
spatial slices, as illustrated in Fig. [T] (see [Zj for details). The integer- valued proper 
time t can be extended in a natural way to the interiors of the four-simplices, 
leading to a global foliation of any causal dynamically triangulated spacetime 
into piecewise flat (generalized) triangulations for any constant real value of t ^2] • 
Inside each building block this time coincides with the proper time of Minkowski 
space. Moreover, it can be seen that in the piecewise linear geometries the mid- 
points of all spatial tetrahedra at constant time t are separated a flxed time-like 
geodesic distance 1 (in units of at) from the neighbouring hypersurfaces at t — 1 
and t+1. It is in this sense that the "link distance" t, i.e. counting future-oriented 
time-like links between spatial slices is a discretized analogue of their proper-time 
distance. 

Let us furthermore assume that the two possible link lengths are related by 

a^ = —aa^. (5) 

All choices a > correspond to Lorentzian and all choices a < —7/12 to Eu- 
clidean signature, and a Euclideanization of geometry is obtained by a suitable 
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analytic continuation in a (see jfT for a detailed discussion of this "Wick rotation" 
where one finds S's(— a) = iSi^{a) for a > 7/12). 

Setting a = —1 leads to a particularly simple expression for the (Euclidean) 
Einstein-Hilbert action of a given triangulation T (since all four-simplices are 
then identical geometrically), namely, 

SEiT) = -koNoiT) + hN,{T), (6) 

with Ni{T) denoting the number of i-dimensional simplices in T. In (jHl), fco is 
proportional to the inverse (bare) gravitational coupling constant, fco ~ I/Cat, 
while k^ is a linear combination of the cosmological and inverse gravitational 
coupling constants. The action (jHl) is calculated from Regge's prescription for 
piecewise linear geometries. If we take a ^ —1 the Euclidean four-simplices of 
type (1,4) and type (2,3) will be different and appear with different weights in 
the Einstein-Hilbert action [7j. For our present purposes it is convenient to use 
the equivalent parametrization 

SEiT) = -koNo{T) + k^N,{T) + A{2Nu{T) + N2s{T)), (7) 

where Ni4{T) and iV23(T) denote the combined numbers in T of four-simplices 
of types (1, 4) and (4, 1), and of types (2, 3) and (3, 2), respectively. The explicit 
map between the parameter A in eq. ((Tj) and a can be readily worked out [Hj . 
For the simulations reported here we have used A in the range 0.4-0.6. 

The (Euclidean) discretized analogue of the continuum proper-time propaga- 
tor (JH) is defined by 

G.„,.„A(T(^no),r(3HT),T) =J2^ ^''^'-''^- (8) 

Terr ^^ 

where the summation is over the set Tt of all four- dimensional triangulations of 
topology S^ X [0, 1] (which we in the following always choose to be S'^) and T 
proper-time steps, whose spatial boundary geometries at proper times and T 
are T^^^(O) and T^^\T). The order of the automorphism group of the graph T 
is denoted by Cr- The propagator can be related to the quantum Hamiltonian 
conjugate to t, and in turn to the transfer matrix of the (Euclidean) statistical 
theory ?]. 

It is important to emphasize again that we rotate each configuration to a 
Euclidean "spacetime" simply in order to perform the summation in the path 
integral, and that this is made possible by the piecewise linear structure of our 
geometry and the existence of a proper-time foliation. Viewed from an inherently 
Euclidean perspective there would be no motivation to restrict the sum over 
geometries to "causal" geometries of the kind constructed above. We also want 



to stress that the use of piecewise hnear geometries has allowed us to write down 
a (regularized) version of (0)) using only geometries, not metrics (which are of 
course not diffeomorphism- invariant), and finally that the use of building blocks 
has enabled the introduction of a diffeomorphism-invariant cut-off (the lattice 
link length a). 

3 Numerical analysis of the model 

While it may be difficult to find an explicit analytic expression for the full prop- 
agator (jH)) of the four-dimensional theory, Monte Carlo simulations are readily 
available for its analysis, employing standard techniques from Euclidean dynam- 
ically triangulated quantum gravity jTHj- Ideally one would like to keep the 
renormalized ^ cosmological constant A fixed in the simulation, in which case 
the presence of the cosmological term A J ^/g in the action would imply that the 
four- volume V4 fiuctuated around (V4) ~ A~^. However, for simulation-technical 
reasons one fixes instead the number A^4 of four-simplices (or^ the four-volume 
V4) from the outset, working effectively with a cosmological constant A ~ V^~^- 

3.1 The global dimension of spacetime 

A "snapshot" , by which we mean the distribution of three-volumes as a function 
of the proper time < t < T for a spacetime configuration randomly picked 
from the Monte Carlo-generated geometric ensemble, is shown in Fig. [21 One 
observes a "stalk" of essentially no spatial extension (with spatial volumes close 
to the minimal triangulation of S^ consisting of five tetrahedra) expanding into 
a universe of genuine "macroscopic" spatial volumes, which after a certain time 
T <T contracts again to a state of minimal spatial extension. As we emphasized 
earlier, a single such configuration is unphysical, and therefore not observable. 
However, a more systematic analysis reveals that fiuctuations around an overall 
"shape" similar to the one of Fig.|21are relatively small, suggesting the existence of 
a background geometry with relatively small quantum fiuctuations superimposed. 
This is precisely the scenario advocated in Sec. 1 and is rather remarkable, given 
that our formalism is background- independent. Our first major goal is to verify 
quantitatively that we are indeed dealing with an approximate four-dimensional 
background geometry ^^ [TTj, and secondly to determine the effective action 
responsible for the observed large-scale features of this background geometry 

UHllIll . 

^For the relation between the bare (dimensionless) cosmological constant k^ and the renor- 
malized cosmological constant A see ^. 

"^For fixed a (or A) one has {N14) ex (A^23) oc {N4). V4 is given as (see for details): 
Vi = ai{NiW&a + 3 + iV23\/12a + 7). We set a^ = 1. 




Figure 2: Snapshot of a "typical universe" consisting of approximately 91000 four- 
simplices as it appears in the Monte Carlo simulations at a given "computer time" . We 
plot the three-volume at each integer step in proper time, for a total time extent of 
T = 40, in units where a^ = 1. 

Important information is contained in how the expectation values of the vol- 
ume Vs of spatial slices and the total time extent r (the proper-time interval 
during which the spatial volumes V3 ^ 1) of the observed universe behave as 
the total spacetime volume V4 is varied. We find that to good approximation 
the spatially extended parts of the spacetimes for various four- volumes V4 can be 
mapped onto each other by rescaling the spatial volumes and the proper times 
according to 

V, - V^/V!'\ r ^ r/Vl'\ (9) 

To quantify this we studied the so-called volume-volume correlator 



(^3(0)1-3(5)) = ^X](l^3(j)^3(j+ 5)) 
^ 3=1 



(10) 



for pairs of spatial slices an integer proper-time distance 5 apart. Fig. El shows 
the volume-volume correlator for five different spacetime volumes V4, using the 
rescaling 0*^, and exhibiting that it is almost perfect. An error estimate yields 
(i = 4 ± 0.2 for the large-scale dimension of the universe [Tij . 

®In (fTn|l we use discrete units sucli tliat successive spatial slices are separated by 1. For con- 
venience we periodically identify T^^^ (T) = T^^^ (0) and sum over all possible three-geometries 
T*^'^^(0), rather than working with fixed boundary conditions. In this way (|10|l becomes a con- 
venient translation-invariant measure of the spatial and temporal extensions of the universe 
(see [H] for a detailed discussion). 
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Figure 3: The scaling of the volume-volume correlator, as function of the rescaled 
time variable x = 5/{NiY/^. Data points come from system sizes A'4 = 22500, 
45000, 91000, 181000 and 362000 at ko = 2.2, A = 0.6 and r = 80. 



Another way of obtaining an effective dimension of the nonperturbative ground 
state, its so-called spectral dimension Ds, comes from studying a diffusion process 
on the underlying geometric ensemble. On a c?- dimensional manifold with a fixed, 
smooth Riemannian metric QabiOy ^^^ diffusion equation has the form 

^ir,(e,eo;^) = A,ir,(e,eo;^), (11) 

where cr is a fictitious diffusion time, A^ the Laplace operator of the metric gabiO 
and Kg{^,^Q;a) the probability density of diffusion from point ^0 to point ^ in 
diffusion time a. We will consider diffusion processes which initially are peaked 
at some point ^o; so that 

^.(e,eo;a=o) = -=i=5^(^-^o). (12) 

Vdet^(0 
For the special case of a fiat Euclidean metric, we have 

p-rf?(C,Co)/4(7 



^,(e,eo;^) 



(47r(T)'^/2 



9abi0=^ab, 



(13) 
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where dg denotes the distance function associated with the metric g. 

A quantity which is easier to measure in numerical simulations is the average 
return probability Pg{cr), defined by 

Pgia) -.= ^1 rf'^e V^d^MO ^.(e, e; ^), (14) 

where V is the spacetime volume V = J d'^^ y^det g{^). For an infinite fiat space, 
we have Pgifj) = l/(47rcr)'^/^ and thus can extract the dimension d by taking the 
logarithmic derivative 

_,dlogPgial^^^ (15) 

a logo- 
independent of a. For nonfiat spaces and/or finite volume V, one can still use 
eq. p5|l to extract the dimension, but there will be correction terms (see jHj for 
a detailed discussion). 

In applying this set-up to four- dimensional quantum gravity in a path integral 
formulation, we are interested in measuring the expectation value of the average 
return probability -Pg(cr). Since -Pg(cr) defined according to (jl4|) is invariant under 
reparametrizations, it makes sense to take its quantum average over all geometries 
of a given spacetime volume V4, 

PvA^) = ^^ hlgab] e-^-^^^-^6{ f d'x,/d^g -V,) Pg{a), (16) 

where Ze{V/i) is the quantum gravity partition function for spacetimes with con- 
stant four- volume V4. 

Our next task is to define a diffusion process on the class of metric spaces under 
consideration, the piecewise linear structures defined by the causal triangulations 
T. We start from an initial probability distribution 

Kr{i,io\(y = Q) = 5i^io, (17) 

which vanishes everywhere except at a randomly chosen (4,l)-simplex io; and 
define the diffusion process by the evolution rule 

KrU, iQ]a + l) = -Y^ Kr{k, i^- a), (18) 

where the diffusion time a now advances in discrete integer steps. These equations 
are the simplicial analogues of (fT^ and (fTT|) . k —>■ j denoting the five nearest 
neighbours of the four-simplex j. In this process, the total probability 

Y,Kr{j,to;(T) = l (19) 
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Figure 4: The spectral dimension Ds of the universe as function of the diffusion 
time 0", measured for kq = 2.2, A = 0.6 and t = 80, and a spacetime volume 
A'4 = 181k. The averaged measurements lie along the central curve, together with 
a superimposed best fit Ds{cf) = 4.02 — 119/(54 + 0") (thin black curve). The two 
outer curves represent error bars. 

is conserved. The probability to return to the simplex zq is then defined as 
-Pr(^o; (^) = Kr{io, io', o") and its quantum average as 






'N4 



1 



(20) 



ioeT]V4 



where T^^ denotes a triangulation with N4 four-simplices, and Se(Xn4) and 
Ze{N4) are the obvious simplicial analogues of the continuum quantities in eq. 

(nni). 

We can extract the value of the spectral dimension Ds by measuring the 
logarithmic derivative as in (jl5|) above, that is. 



Dsia) 



dloga 



(21) 



as long as the diffusion time is not much larger than N/ '^ . The outcome of the 
measurements is presented in Fig. EJ with error bars included. (The two outer 
curves represent the envelopes to the tops and bottoms of the error bars.) The 
error grows linearly with a, due to the presence of the logo" in (j?H) . 
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The remarkable feature of the curve Ds((t) is its slow approach to the asymp- 
totic value of Ds{cr) for large a. The new phenomenon we observe here is a scale 
dependence of the spectral dimension, which has emerged dynamically [T7l[T3j. 

As explained in [T7], the best three-parameter fit which asymptotically ap- 
proaches a constant is of the form 

Ds{(T) = a ^ = 4.02--^^. (22) 

The constants a, b and c have been determined by using the data range a G 
[40, 400] and the curve shape agrees well with the measurements, as can be seen 
from Fig. EJ Integrating ()22j) we obtain 

Pi^) - -TFTjA-n-^.^ (23) 



a"/2(l + c/a)^/' 



from which we deduce the limiting cases 

^-«/2 fQj, large (J, 

P{a) ~ <; (24) 

Again we conclude that within measuring accuracy the large-scale dimension of 
spacetime in our model is four. We also note that the short- distance spectral 
dimension seems to he approximately Dg = 2, signalling a highly non-classical 
behaviour. 



3.2 The effective action 

Our next goal will be to understand the precise analytical form of the volume- 
volume correlator l\l(j\} . To this end, let us consider the distribution of differences 
in the spatial volumes V3 of successive spatial slices at proper times t and t + 6, 
where S is infinitesimal, i.e. 5 = 1 in lattice proper-time units. We have measured 
the probability distribution Py^ (z) of the variable 

for different values of V3. As shown in Fig. they fall on a common curve. ^ 
Furthermore, the distribution Pvii^) is fitted very well by a Gaussian e~^^ , with 
a constant c independent of V3. From estimating the entropy of spatial geome- 



^Again we have applied finite-size scaling techniques, starting out with an arbitrary power 
Kj" in the denominator in H25|) . and then determining a = 1/2 from the principle of maximal 
overlap of the distributions for various Va's. 
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Figure 5: Distribution Py^i^z) of volume differences of adjacent spatial slices, for three- 
volumes Vz = 10.000, 20.000, 40.000 and 80.000 tetrahedra. 



tries, that is, the number of such configurations, one would expect corrections of 
the form V3", with < a < 1, to the exponent cz^ in the distribution Py^{z). 
Unfortunately it is impossible to measure these corrections directly in a reliable 
way. We therefore make a general ansatz for the probability distribution for large 
V^{t) as 



exp 



^3(t) 



dt 



-CaK^W 



(26) 



where < a < 1, and C\ and Ci are positive constants. 

In this manner, we are led by "observation" to the effective action 



S. 



V4 



dt 



Cl 



^3(t) 



dV^3(t) 

dt 



c-,V^{t)-\V^{t) 



(27) 



valid for large three- volume V3(t), where A is a Lagrange multiplier to be deter- 
mined such that 



dtV^{t) = V^. 



(21 



From general scaling of the above action it is clear that the only chance to obtain 
the observed scaling law, expressed in terms of the variable t/V^ ? is by setting 

a = 1/3. In addition, to reproduce the observed stalk for large times t the 

1/3 
function V^ has to be replaced by a function of V3 whose derivative at goes 
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like V^, z/ > 0, for reasons that will become clear below. A simple modification, 
which keeps the large- V3 behaviour intact, is given by 

V^V3 ^ (1 ^ y^y/s _ 1^ (29) 

but the detailed form is not important. If we now introduce the (non-negative) 
scale factor a(t) by 

Vsit) = a\t), (30) 

we can (after suitable rescaling of t and a(t)) write the effective action as 

Sy^ = ^[^^ U^) (^) ' + <^) - ^«' W) ' (31) 

with the understanding that the linear term should be replaced using (jHUj) and 
()29|1 for small a{t). We emphasize again that we have been led to (j!?T|) entirely 
by "observation" and that one can view the small-a(t) behaviour implied by (j29|) 
as a result of quantum fluctuations. 

3.3 Minisuperspace 

Let us now consider the simplest minisuperspace model for a closed universe 
in quantum cosmology, as for instance used by Hartle and Hawking in their 
semiclassical evaluation of the wave function of the universe |T9j. In Euclidean 
signature and proper-time coordinates, the metrics are of the form 

ds^ = df + a^{t)dnl, (32) 

where the scale factor a{t) is the only dynamical variable and dQ^ denotes the 
metric on the three-sphere. The corresponding Einstein-Hilbert action is 

S'' = ^J d^ {-<t) (^) ' - <t) + Aa^(^)] • (33) 

If no four- volume constraint is imposed, A is the cosmological constant. If the 
four-volume is fixed to V4, such that the discussion parallels the computer simu- 
lations reported above, A should be viewed as a Lagrange multiplier enforcing a 
given size of the universe. In the latter case we obtain the same effective action 
as that extracted from the Monte Carlo simulations in (j!?T|l . up to an overall sign, 
due to the infamous conformal divergence of the classical Einstein action evident 
in ()33p. From the point of view of the classical equations of motion this overall 
sign plays of course no role. Let us compare the two potentials relevant for the 
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calculation of semiclassical Euclidean solutions associated with the actions ()33|) 
and dSD). The "potential" i° is 



V{a) 



-a + \a 



(34) 



and is shown in Fig. El without and with sniall-a modification, for the stan- 
dard minisuperspace model and our effective model, respectively. The quantum- 





Figure 6: The potential V{a) of (|34|) underlying the standard minisuperspace dynamics 
(left) and the analogous potential in the effective action obtained from the full quantum 
gravity model, with small-a modification due to quantum fluctuations (right). 

induced difference for small a is important since the action ()3H) admits a clas- 
sically stable solution a{t) = which explains the "stalk" observed in the com- 
puter simulations (see Fig. |2)). Moreover, it is appropriate to speak of a Euclidean 
"bounce" because a = is a local maximum. If one therefore naively turns the 
potential upside down when rotating back to Lorentzian signature, the metastable 
state a{t) = can tunnel to a state where a{t) ~ I4 , with a probability ampli- 
tude per unit time which is (the exponential of) the Euclidean action. 

In order to understand how well the semiclassical action (j31|) can reproduce 
the Monte Carlo data, that is, the correlator (fTI]|) of Fig. |2l we have solved for the 
semiclassical bounce using (jSU, and presented the result as the black curve in Fig. 
El The agreement with the real data generated by the Monte Carlo simulations 
is clearly perfect. 

The picture emerging from the above for the effective dynamics of the scale 
factor resembles that of a universe created by tunneling from nothing (see, for 

^°To obtain a standard potential - without changing "time" - one should first transform to a 
variable x = a^ for which the kinetic term in the actions assumes the standard quadratic form. 



It is the resulting potential V{x) 
for small x such that V^'(O) = 0. 



= _.^2/3 



+ Ax^ which in the case of (|^ should be modified 
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example, [2011^1221); although the presence of a preferred notion of time makes 
our situation closer to conventional quantum mechanics. In the set-up analyzed 
here, there is apparently a state of vanishing spatial extension which can "tunnel" 
to a universe of finite linear extension of order a ~ V^ . Adopting such a 
tunneling interpretation, the action of the bounce is 



V, 



1/2 



St - 7^' (35) 

and the associated probability per unit proper time for the tunneling given by 



-s. 



<!# 



P{Vi) ~ e"''^4 . (36) 



4 Discussion 



Causal dynamical triangulations (CDT) provide a regularized model of quantum 
gravity, which uses a class of piecewise linear geometries of Lorentzian signature 
(made from fiat triangular building blocks) to define the regularized sum over 
geometries. The model is background-independent and has a diffeomorphism- 
invariant cut-off. For certain values of the bare gravitational and cosmological 
coupling constants we have found evidence that a continuum limit exists. The 
limit has been analyzed by rotating the sum over geometries to Euclidean sig- 
nature, made possible by our use of piecewise linear geometries. The geometries 
included in the sum thus originate from Lorentzian-signature spacetimes, a class 
different from (and smaller than) the class of geometries one would naturally in- 
clude in a "native" Euclidean path integral. We have concentrated on computing 
a particular diffeomorphism-invariant quantity, the proper-time propagator, rep- 
resenting the sum over all geometries whose space-like boundaries are separated 
by a geodesic distance T. The sum over such geometries allows a simple and 
transparent implementation in terms of the above-mentioned building blocks. 

In the Euclidean sector of the model, which can be probed by computer sim- 
ulations we observe a four-dimensional macroscopic universe that can be viewed 
as a "bounce". When we integrate out (a/f:er having constructed the full path 
integral) all geometric degrees of freedom except for the global scale factor, the 
large-scale structure of the universe (the bounce) is described by the classical 
general-relativistic solution for a homogenous, isotropic universe with a cosmo- 
logical constant on which (small) quantum fluctuations are superimposed. We 
find this result remarkable in view of the difficulties - prior to the introduction 
of causal dynamical triangulations - to obtain dynamically any kind of "quan- 
tum geometry" resembling a four- dimensional universe. In our construction, the 
restrictions imposed by causality before rotating to Euclidean signature clearly 
have played a pivotal role. 



A number of issues are being addressed currently to obtain a more complete 
understanding of the physical and geometric properties of the quantum gravity 
theory generated by CDT, and to verify in more detail that its classical limit is 
well defined. Among them are: 

(i) A better understanding of the renormalization of the bare coupling constants 
in the continuum limit, with the currently favoured scenario being that of asymp- 
totic safety |2S1- There are very encouraging agreements between the results of 
CDT and those of a Euclidean renormalization group approach |24] (see [25j for 
older, related work). In particular, both approaches obtain a scale-dependent 
spectral dimension which varies between four on large and two on short scales, 
(ii) An identification and measurement of the "transverse" gravitational degrees 
of freedom, to complement the information extracted so far for the scale factor 
only. For background-independent and coordinate-free formulations like CDT we 
still lack a simple and robust prescription for how to extract information about 
the transverse degrees of freedom, a quantity analogous to the Wilson loop in 
non-abelian gauge theories. 

(iii) The inclusion of matter fields in the computer simulations. Of particular 
interest would be a scalar field, playing the role of an inflaton field. While it is 
straightforward to include a scalar field in the formalism, it is less obvious which 
observables one should measure, being confined to the Euclidean sector of the 
theory. Based on a well-defined CDT model for the nonperturbative quantum 
excitations of geometry and matter, moving the discussion of quantum cosmol- 
ogy and various types of inflation from handwaving arguments into the realm of 
quantitative analysis would be highly desirable and quite possibly already within 
reach. 
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